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Abstract 



A scenario for the generation of the baryon asymmetry in the early Uni- 
verse is proposed in which cosmic string loops, predicted by theories where 
the baryon and/or lepton numbers are gauged symmetries, collapse during 
the friction dominated period of string evolution. This provides a mecha- 
nism for the departure from thermal equilibrium necessary to have a non- 
vanishing baryon asymmetry. Examples of models are given where this idea 
can be implemented. In particular, the model with the gauge symmetry 
SU(3) c <g)SU(2) L <g)U(l)Y ®U(i) B ®U(i) L has the interesting feature where 
sphaleron processes do not violate the baryon and lepton numbers so that no 
wash out of any initial baryon asymmetry occurs at the electroweak scale. 
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1. Introduction 

The basic requirements for the generation of the baryon asymmetry of the 
Universe (BAU) |I| are: (1) baryon number violation, (2) C and CP violation 
and (3) a departure from thermal equilibrium. In the past, models of Grand 
Unified Theories (GUTs) in conjunction with the standard hot big model of 
cosmology were found to have the above necessary ingredients to explain the 
observed matter-antimatter asymmetry of the Universe . However typical 
GUTs require relatively high symmetry breaking scales (usually not too far 
away from the Planck scale) making direct testing of such theories inaccessi- 
ble to currently known experimental techniques. It would be interesting to 
find particle physics models which can realise the baryogenesis scenario and 
can at the same time be tested experimentally in the near future. Recent at- 
tempts to do this include baryogenesis at the electroweak scale using just the 
standard model (SM) of particle physics ||. In this paper we will explore an 
alternative mechanism for generating the BAU and give examples of particle 
physics models which can implement such a scenario. 

This mechanism of generating the BAU is based on the collapse of cosmic 
string loops. A variant of this mechanism has been studied in Ref. ||. This 
will be described in section 2. Then in the following sections we will give 
examples of particle physics models which can realise this mechanism. The 
first example will be a model with gauged baryon and lepton numbers. The 
crucial aspect of this model is that since the baryon and lepton symmetries 
are anomaly free then the sphaleron processes will not violate these symme- 
tries and hence any initial asymmetry produced will not be washed out. This 
is in contrast with many of the GUT scenarios for the generation of the BAU 
where the wash out of the initial asymmetry is a common problem. The 
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second example we will consider is a model where the (B — L) symmetry is 
gauged. The right-handed neutrinos in this model are responsible for estab- 
lishing an initial lepton number asymmetry. The lepton number asymmetry 
is then converted to a baryon number (£>) asymmetry via the sphaleron. 
In both of these models the ultimate source of the asymmetry comes from 
particle production as a result of the collapse of cosmic string loops. 

We conclude this section by commenting on the motivation for the par- 
ticular particle physics models we have chosen to discuss. We feel that these 
models are interesting simply from the particle physics point of view, let alone 
their implications for cosmology. Firstly, these models deal with physics be- 
yond the SM at scales well below the GUT scale. They are simple in the 
sense that they rely only on the concepts central to the theoretical frame- 
work of the SM. They are (i) gauge symmetry and (ii) spontaneous symmetry 
breaking. The success enjoyed by the SM in describing nature up to energies 
hitherto explored by experiment has promoted these two concepts to a status 
of fundamental importance. So much so, that in the case of symmetries, any 
symmetry which is not gauged is viewed with suspicion by some. For exam- 
ple, even though global symmetries (continuous or discrete) appear in many 
particle physics models, including the SM, they are considered to be less than 
fundamental or at best "accidental". Since the origin of symmetries is still 
an open question, such a theoretical prejudice might prove to be wrong in the 
future. Nevertheless, in this paper we will pursue this theoretical prejudice, 
perhaps to an extreme, by gauging all or some combination of the known 
classical global symmetries of the SM. 
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2. Collapsing string loops mechanism 

In this section we will describe a mechanism by which the collapse of 
cosmic string loops can provide a source of baryon or lepton number. The 
analysis will have some similarities to that given in Ref. 0] but because of 
the different scales involved there will be some important differences. 

Consider the Abelian Higgs model (which is the simplest example of a 
model in which cosmic strings can form as a result spontaneous symmetry 
breaking ||) where a complex Higgs field develops a vacuum expectation 
value (VEV) (H) = ae 10 . The arbitrary phase, 9 = 9(x), can vary in different 
regions of space. For H to be single valued 9 must change by an integer 
multiple of In around a closed loop. When the loop is shrunk to a point 
9 becomes undefined so that there exists a point where (H) = 0, i.e., thin 
tubes of false vacuum get trapped somewhere inside the loop. Such tubes of 
false vacuum or "strings" must either be closed or infinite in length so that 
the closed loop cannot be contracted to a point without encountering the 
tubes of false vacuum. Hence these strings are topologically stable. 

Once the cosmic strings have been produced their evolution in time is 
demarcated by two periods. Initially the strings will be moving in a dense 
medium so that their motion is strongly damped and this is usually referred 
to as the friction dominated period. After this time the friction becomes 
negligible and the strings can move at relativistic speeds. The evolution of 
strings is then mainly determined by the expansion of the Universe, grav- 
itational radiation of string loops and the intercommuting of intersecting 
strings. The friction dominated period lasts from the time of the phase tran- 
sition, t c , until the time t* ~ (G/i) -1 ^, where G is Newton's constant and 
fi ~ cr 2 is the mass per unit length of the string. In terms of temperatures, 
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T* ~ jfi-T c , where M p is the Planck mass and T c ~ a is the critical tempera- 
ture of string formation. For instance, if the phase transition takes place at 
scales less than a ~ 10 10 GeV then the friction dominated period will not 
end before the temperature drops below about 10 GeV. Since in this paper 
we are interested in relatively low symmetry breaking scales, only the friction 
dominated period of string evolution will be relevant. The situation depicted 
here is therefore different from that described in Ref. [§J where the breaking 
scale is very close to the GUT scale and the friction dominated period lasts 
for a much shorter time. In such the typical scale is the coherence 

length for a long string network and the fact that its expansion is catching 
up with the Hubble radius must be used. 

Also note that since the gravitational effects of cosmic strings depend 
on the dimensionless quantity G/i, and a -C M p , then we expect that the 
usual cosmological effects associated with strings (e.g. seeding for structure 
formation) will be negligible. 

The basic idea for a nonzero B and/or L produced from strings is that 
when strings collapse to a size comparable to its width the microphysical 
forces of the string become important. As a result the energy of the string will 
be released in the form of particle production. We will assume for simplicity 
that all the energy released from the collapsed string goes into the quanta 
of the massive Higgs field. The energy of Nq quanta of mass ~ X^a 
produced from the collapse of a single string loop is NgX^a, where A is the 
quartic self-coupling of the Higgs field. The energy of the string loop is given 
by nf3R where f3 ~ 2ti and R is the loop radius. Equating the two expressions 
for the energy for R ~ w results in Nq ~ /3A _1 , where the thickness of the 
string w ~ X^^a^ 1 was used. 

We will now proceed to estimate the number of loops contributing to 
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particle production. The size distribution of loops, for loops of size I, is 

dn ~ £- 4+a r a dl, (1) 

where £ c ~ (Act) -1 is the correlation length of the Higgs field and a was found 
from Monte Carlo simulations of string systems to take on values from 2 to 
2.5 0. Integrating Eq. ([[J) gives 

n{t w ) ~ {/(t c ) (1 ~ a) - , (2) 

where t w is the time when l{t w ) = w and l(t c ) is expected to be the largest 
sized loop that can collapse to l(t w ) in time. 

To find the explicit time dependence of the loop size, consider the motion 
of the string in the friction dominated epoch. The typical damping time on 
the motion of the string is J7] 

AP/l a a 2 

where the change in momentum per unit length of the string is AP/l ~ fiv 
and the force per unit length on the string is F/l ~ nv . v is the speed of 
the string relative to the surrounding medium and n = k n T 3 is the number 
density with 

*» = ^ £ *f§) 3 + ! £ S <@ 3 . (4) 
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i=bosons i=fermions 

where ((x) is the Riemann zeta function and gi is the number of degrees of 
freedom of the particles in the dense medium. Therefore the characteristic 
loop size at a time t is 



where the time-temperature relation of t = nM p T~ 2 for an expanding Uni- 
verse in the radiation dominated epoch was used. The rate at which the 
loops collapse can be estimated by assuming 

dl dtd 
dt dt 

This then gives 

'2 
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with l(t w ) = w. Substituting this into Eq. (0) results in 
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This is the number of loops per unit volume that will collapse to size w. 

However, not all of these loops will contribute to generating a nonzero 
B or L since the decays of the Higgs particles will be compensated by their 
inverse decays at temperatures above the freeze-out temperature, Tp, (cor- 
responding to a time, tp) which is of the order of the Higgs boson mass. The 
relevant number density of loops is then given by 



n 



tF 



a(t v 



where 
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dn(t, u 



dt 



(9) 
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and a(t) oc oc T~ l is the usual cosmic scale factor in the radiation dom- 
inated epoch of the Universe. The limits of integration correspond to their 
respective temperatures, i.e., 

_i i 

train < > T m i n — kn 3 A 6 (J 



t F 



A 2 a, 



(11) 
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T min is obtained by requiring l(t) > w (see Eq. (||)). To simplify the evalua- 
tion of n, notice that the ratio 



a{t T 



1 3 / rp \ 3 

-LF 



a(tj 



T 
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KX. (12) 



This is expected to be of order one for "typical" values of the parameters 
(e.g. k n ~ 10 and A ~ 0.1). Therefore the scale factor ratio shouldn't vary 
too rapidly for t F < t < t min and hence will be neglected. Then Eq. (f]) 
becomes 



Tl — _ 

(1 - a) [ V &„A, 

So the number density of heavy particles produced after tp is given by Nqu 
and this leads to a baryon number (or lepton number depending on the model 
under consideration) given by 

^ / t \ Nan 
B ( or L) ~ — =-e 



27r 2 p* (1 - at) [ V £; n A, 

where e is the particle physics model dependent CP-violation factor and 



s ~ j^^ 2 g* Aacr 3 is the entropy density evaluated at T F . 
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3. Gauged baryon and lepton numbers 

From an experimental point of view, baryon and lepton number appear to 
be good symmetries, at least at the classical level. An interesting possibility 
for physics beyond the SM is to gauge the baryon and/or lepton numbers. 
However, no long distance effects associated with baryon or lepton number 
have been observed. So if they are gauge symmetries, then they must be 
spontaneously broken. 
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To construct a model where B and/or L are gauged, the fermion sector 
of the SM needs to be extended so that these symmetries are anomaly-free. 
We know of two ways to do this. One way is to add a mirror set of fermions 
which have the same quantum numbers as the SM fermions but with opposite 
chiralities. In this case the anomalies cancel trivially. The other way is to 
add "exotic" generations to the SM fermions. The multiplet structure of an 
exotic generation can be defined as: 

leptons : N R {Y) © (JV - 1) L (Y - 1) © (JV - 1)' L (Y + 1) © (JV - 2) R (Y); 
quarks: N R (V) © (JV - 1) L (~Y - l) 

®(N - i)' L (-hy + 1) © (jv - 2) R (-~y) , (15) 

where N is the JV-dimensional representation of SU(2)l and Y is the hyper- 
charge. One of the simplest anomaly-free structures of the gauge group 

G BL = SU(3) C ® SU(2) L ® U(1) Y ® U(1) B ® U(1) L (16) 

corresponds to four ordinary doublet-type and two exotic triplet-type gener- 
ations. Note that it is also possible to use the same methodology to construct 
models where either B or L are anomaly-free only. For further details see 
Ref. gj. 

Having noted that it is possible to gauge B and/or L, we now consider a 
model with the gauge group given by Eq. ([L6D where we will apply the re- 
sults of section 2 to generate the BAU. The fermion sector consists of the SM 
fermions together with a set of either mirror fermions or exotic generations 
to cancel the potential B and L anomalies. The usual electroweak symmetry 
breaking and fermion mass generation are facilitated by the SM Higgs dou- 
blet. In addition, two complex SM Higgs singlets, Hi and H 2 , are introduced 
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to spontaneously break U(1)b and U(1)l when they acquire nonzero VEVs, 
(ae ld ) and (rje' 1 ^), respectively. The Higgs fields transform under the gauge 
group G B l as 

fT 1 ~(l,l,0)(2,0) and H 2 ~ (1, 1, 0)(0, 2), (17) 

where the first set of parentheses give the SM quantum numbers and the sec- 
ond set gives the B and L numbers. From phenomenological considerations, 
it is expected that the scales of B and L breaking lie above the electroweak 
scale. Since B and L are spontaneously broken the formation of cosmic 
strings will result and particle production from the collapse of these string 
loops will follow the scenario described in section 2. In this section we will 
only be interested in what happens in the quark sector for the generation of 
the BAU and thus we will not consider the lepton sector any further. So far 
a nonzero VEV for Hi only breaks the gauged B symmetry, leaving behind 
a residual global B symmetry. To also break the global symmetry so that 
there is baryon number violation in the quark sector, the following vectorlike 
pair of fermions are introduced. Their transformation properties under Gbl 
are defined by their Yukawa interactions with Hi. 

C$ uk = XlHiQ l Q c l + X R HiQ R Q c R + 2DQ L Q R + H.c, (18) 

so that 

Q L ~(8,1,0)(1,0) and Q R ~ (8, 1,0)(1,0). (19) 

Finally, we need to introduce one more scalar field into the model. This is 
a charged colored scalar field, A, which does not get a nonzero VEV. Its 
purpose is to transfer the B asymmetry produced in the Ql, R sector to the 
ordinary quark sector. The transformation properties of A are determined 
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by the interaction term, 

Cyuk = X A Q L Ad R + H.c, (20) 

so that A ~ (3, 1, §)(§, 0), where cLr is the usual right-handed (i-type quark 
field. 

We will now proceed to determine the BAU from the collapse of string 
loops by using the results of section 2. First we will outline the scenario 
of section 2 as applied to this model. When the Universe cools down to a 
temperature, T c ~ a, the U(1)b symmetry is spontaneously broken and the 
formation of cosmic strings is initiated. The string loops collapse rapidly in 
the friction dominated epoch of the evolution of the string network. The 
production of Higgs scalars, h (h is a real scalar field left over from the 
Higgs mechanism with its mass given by ~ A^cr), is assumed to result 
from this collapse. This provides a source of these Higgs scalars which are 
at this point overabundant and hence out of equilibrium. The scalar, h, 
then decays rapidly into the Q-type particles since its decay width, of the 
order of M^, is much larger than the Universe expansion rate, H ~ T 2 /M p , 
at the temperatures under consideration and thereby an excess of B can 
be generated. The Q's in turn get converted to d-type quarks via their 



interactions with the A scalars. B can be calculated from Eq. (14) once e is 
determined. Hence the remainder of this section is used to show how e can 
be calculated. 

To determine e we need to write down the interactions of h with the Q's. 
First write Ql,r in terms of their mass eigenstates. When H\ develops a 
nonzero VEV the interaction terms of Eq. fllTf ) generates the following mass 
terms: 

C mass = VrM-Vl + H.c, (21) 
10 



where 

•*-(£)• -(£)■ (22 > 

Note that A^ can be made real by a phase redefinition of the fields. The 
mass matrix, A4, can be diagonalized so that 

C ma .. = *' R TW , L + &c., (23) 

where V = Diag(Mi, M 2 ) and the mass eigenstates are given by 

*i=(^)=^. **=(S)= yt * fl (24) 

with 

tt — ( cos $ si n ^ \ j T/ _ f cos6> sin^ \ , . 

\ -e i5 sm6 e i5 cos6 ) = { -e" w sin0 e^cosfl J ' ( j 

The interactions of ft from Eq. (|1|) can now be written as 

dm = \iK\ l Cil + \2K c m C2R + XwKmtiL + H.c, (26) 

where 

Ai = Xl cos 2 6 + \Re i2S sin 2 6, 
A 2 = Ai sin 2 6 + A^e - * 25 cos 2 6, 

A 12 = sin2^(A L -A fi e i25 ). (27) 

It is important to note that Ai j2 and Ai 2 cannot be made real because the 
phases of the fields are now fixed in the mass eigenstate basis. These complex 
couplings are necessary for the nonzero CP-violation parameter, e, to be 
generated. From Eq. fl26|) , it can be seen that there are three processes 
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which contribute to e. They are (1) h — > (il&r, (2) h — > CilCil an d (3) 
h — ► (2r(2R- The Feynman diagrams for process (1) are shown in Figs. 
(1) and (2). The corresponding diagrams for processes (2) and (3) are very 
similar. For example, for process (1), only the diagram in Fig. 2(b) has a 
complex product of couplings for the interference term between the tree and 
one-loop amplitudes and hence can contribute to e. The relevant amplitudes 
can be written as follows: 

A(h — > (il(2r) = Ai 2 So + A^AiA^ai 

A{h^(iL(2 R ) = A; 2 oo + Ai 2 AIA 2 ai, (28) 

where a and a\ denote the tree and one-loop contributions respectively. 
Then 

\A\ 2 - \A\ 2 = -4 Im (A1A2A12) Im (a*ax) . (29) 
For convenience, define 

A(i — /) = r(< — /) - r(i — /), (30) 

where i and / represent the initial and final states of the decay process. Then 

A(h — > Ci^) , , 4 Im (AiAjAjf) Im (agaQ 3 
r(/i — all) " (|Ax| 2 + |A 2 | 2 + 2|A 12 | 2 ) \atf P 12 1 ' 

Similarly for processes (2) and (3) 

A(/i — > CilCil) _ , 2 Im (AiAjA*!) Im (a* 0l 



r(/i — > all) (|Ax| 2 + |A 2 | 2 + 2|A 12 | 2 ) |a | 2 



P 3 n (32) 



and 



A(fe — > C 2 rC 2 k) „ 2Im(A 1 A 2 ^ 2 )Im(q*q 1 ) 3 
r(^— .all) " (|A 1 | 2 + |A 2 | 2 + 2|A 12 | 2 )|a r 22 1 J 

respectively, where a and a x are the amplitudes a and a\ with the phase 



space factor = Jl — (Mi + Mj) 2 /Ml factored out. Note that a and a± 
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are not the same for each of the above processes since Mi ^ M 2 in gen- 
eral. However in the following, to simplify the calculations, we will assume 
Mi ~ M 2 and much smaller than M^. There is also a set of diagrams where 
the internal Higgs line is replaced by the massive gauge boson line but these 
diagrams do not contribute to e since the product of couplings of the in- 
terference terms are real due to the real gauge couplings. So far, all of the 
S-violation is carried by £1,2 in the exotic sector but this will be transferred 
to the ordinary sector via the interactions of Eq. fl2Cf ), i.e., 



C-Yuk = Aa (cos 6( 1L Ad R + sin 9( c 2R Ad R ) + H.c, (34) 



Therefore e can be written as 

Ajh — > f) vij — » n 

^ f Y(h — >all)r(/ — ^ all) 



B 



f 



+ 



T(h — ► all) r(Ci£ — > all) 

Mh — > c 2fi c 2fi ) r(c 2fl a 1 "^) 

r(/i — ► all) r(C 2 R — ► all) 



1 A(h — > CirCgg) 

2 r(/i — »■ all) 

which simplifies to 



r(CiL — Atd fl ) , T(( 2R — > Atdj 



T(Cil — > all) T(C 2 ii — > all) 



(35) 



e ~ 



1 (Mi — M 2 f Im(AiA^A^) 



(36) 



4tt M 2 |Ai| 2 + |A 2 | 2 + 2|Ai 2 | 2 ' 

where the branching ratios of the processes, Ci£> (zr — > A^dn are unity, the 
factor 167rlm (ag a i) /l°o| 2 is of order one || and Bd = — | for the final state 
baryon number since each process results in the decay product d R . To 
get a rough estimate of the size of the couplings required for the observed 
value of the BAU, let A^ = Ai ~ A 2 ~ A12 and M i 2 ~ O(X^a), so that 
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e ~ \y (167rA). Substituting this in Eq. (|I~4| ) gives 

B ~ 1 ^^{^-( 1 + £) <1 '1* (37) 

For example, by taking A ~ 0.1, a = 2.5, with k n ~ 10 (see Eq. (^)), then 
A^ ~ 10~ 2 for the present baryon asymmetry which lies in the range (4 — 
. Therefore if M 12 ^5 10 3 GeV, then the breaking scale for B 
is expected to be less than about 10 5 GeV. Note that no particular limit must 
be set on the L breaking scale since now sphalerons violate neither B nor L 
and therefore the interplay between the sphalerons and L- violating processes 
is completely harmless for the generated baryon asymmetry. Furthermore, 
the dangerous 5-violating processes like A + A — > d R + d R mediated by 
the ( particles, whose rate is ~ X\T 3 /m^, are out of equilibrium for small 
enough values of Aa- As a result they also do not affect the B produced from 
the rapid decays of ( — > A + d R . 

4. The gauged (B — L) model 

It is known that the SM with right-handed neutrinos has an anomaly-free 
{B — L) symmetry. Therefore (B — L) can be gauged. In this section we will 
consider such a model for generating the BAU using the collapse of string 
loops mechanism of section 2. The gauge group is 

SU{$) c ®SU{2) l ®U{1)y®U(1) B -l (38) 

To break the (B — L) symmetry and to give the right-handed neutrinos a 
Majorana mass, we introduce the scalar field, Xi which is a complex SM 
Higgs singlet with two units of (B — L). As the temperature in the early 
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Universe drops to about T ~ (x) the (B — L) symmetry is spontaneously 
broken and the formation of cosmic strings begins. These strings will evolve 
in the friction dominated period where string loops tend to collapse rapidly. 
The collapse of string loops will result in the production of Higgs particles 
associated with x- These Higgs particles will then in turn decay into right- 
handed neutrinos providing a source of lepton number. By using the results 
of section 2, L will be given by Eq. fll4|). The e parameter can be calculated 
from the decays of the right-handed neutrino, i. e. 



such that 



(39) 



(40) 



where e« is the difference between particle-antiparticle branching ratios given 
by H 

with 



/(*) 



1 - (1 + x)ln 



1 + x 



x 



(42) 



The h's denote the Yukawa couplings between right-handed and left-handed 
neutrinos through the SM Higgs doublet $ and are assumed to be complex 
to give a source of CP violation. 

The lepton number produced by right-handed neutrino decays is con- 
verted to a nonzero baryon number via the (B — L) conserving sphaleron 
processes. Therefore the BAU is given by 



B 



k L 
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^^-'f^l.^'f-l, (43) 



27T 2 <7* (1 — a) { \ knX, 

where k is a numerical factor of 0(1) and can be easily calculated from 
Ref. [ 12 1 . For A ~ 10~ 2 and a = 2.5 gives a baryon asymmetry of about 
B ~ 10 _1 e. For the present baryon asymmetry to be of the order of5xl0~ n , 
e has to be as large as 1CT 10 . 

We also need to check that the initial baryon asymmetry generated by the 
collapsing string loops and the subsequent decays of right-handed neutrinos 
is not erased by a combination of other lepton violating interactions and 
sphaleron processes ||13|| . The lowest dimension L- violating operator in the 
effective low energy theory is given by 

Tfl 

Cal=2 = -^F L F L $$ + H.c, (44) 
ir 

where m Vh is the mass of the left-handed neutrino and v is the electroweak 
breaking scale. The interaction rate of these AL = 2 processes is Tal=2 — 
m1 L T 3 / (7r 3 f 4 ). For the survival of the pre-existing asymmetry we require 
this interaction rate to be less than the expansion rate of the Universe, i.e., 
H ~ T 2 /M p , where T is given by T m i n in Eq. (11). This results in a bound 
on the (B — L) breaking scale, 

(X) > -^n V3 A 1/2 M p , (45) 

where h is the Yukawa coupling for the Dirac mass term and 7 is the one 
for the Majorana mass. So for natural values of the couplings involved, the 
bound on (x) can easily be made as low as 10 9 GeV. 



5. Conclusions 
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In the present paper we have shown that models, in which B and L 
are gauged symmetries, are interesting not only from the particle physics 
perspective but also for cosmology - the generation of the baryon asymmetry 
in the early Universe. Indeed, such models provide a natural implementation 
of the out-of-equilibrium condition through the collapse of string loops formed 
during the phase transition at intermediate scales where B and/or L are 
spontaneously broken. Also, differently from what happens in ref. ||], these 
loops decay during the friction dominated epoch of string evolution. We 
would also like to stress that, in the case in which both B and L are gauged, 
the sphaleron processes do not violate B + L so that there is no wash out 
of any pre-existing asymmetry and no particular limit must be imposed on 
the scale of L breaking (other than that from experiment). Finally, in our 
scenario, the electroweak phase transition need not be first order as required 
in most non-GUT scenarios for the generation of the BAU |J|. 
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Figure 1: Tree level diagram for the decay h — > CilC2R- 
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Figure 2: One-loop diagrams contributing to CP violation in the h — > (il&r decay; 
only diagram (b) gives a nonvanishing contribution. 
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